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A procedure is outlined for the determination of magnetic exchange constants in anisotropic 
perovskite anitferromagnets using powder inelastic neutron scattering. Spin wave densities-of-states 
are measured using time-of- flight inelastic neutron scattering for LaMnC>3 (A-type antiferromagnet), 
LaVO;} (C-type), and LaFeOa (G-type) and compared to Heisenberg model calculations. The 
anisotropy of in-plane (J a b) an d out-of-plane (J c ) exchange constants can be obtained from the 
data. The procedure quickly determines the magnetic exchange interactions without the need for 
single-crystal dispersion measurements and allows for rapid systematic studies of the evolution of 
magnetism in perovskite systems. 

PACS numbers: 75.30.Ds, 75.30.Et, 78.70.Nx 



I. INTRODUCTION 

Many important magnetic materials, such as colossal 
magnetoresistive manganites and high temperature su- 
perconductors, are based on perovskite transition metal 
oxides. The ground states and phases of these materi- 
als are known to depend sensitively on the energy scales 
of magnetic, orbital, vibrational, and electronic degrees- 
of-freedomi^ Various spectroscopic techniques are em- 
ployed to determine these energy scales, and the coupling 
between them, in an effort to understand and control the 
myriad of physical properties of these compounds. 

The magnetic energy scale is set by the exchange en- 
ergy between magnetic ions. It is purely quantum me- 
chanical in origin and can arise from many different pro- 
cesses originating from the exchange of electrons between 
magnetic ions; such as direct exchange, superexchange, 
and double exchange. In insulating transition metal ox- 
ide materials, the superexchange interaction depends on 
the overlap of metal cZ-orbitals on neighboring sites via 
oxygen ligands. The rules for determining the sign (fer- 
romagnetic or antiferromagnetic) and strength of the su- 
perexchange between neighboring ions were established 
early on^ A brief summary of these rules as applied to 
perovskites is as follows; two half-filled (or empty) or- 
bitals in a (180°) bonding configuration have antiferro- 
magnetic (AF) exchange, while a half-filled and an empty 
orbital in a bonding configuration have ferromagnetic (F) 
exchange. 

In cubic perovskites, F exchange is not expected since 
all ionic sites are equivalent. Such is the case for LaFcC>3, 
where each Fe 3+ ion has a half-filled 3d 5 configuration 
and all exchange interactions are AF, leading to a G- 
type magnetic structure (with all neighboring magnetic 
ions aligned antiparallel) . However, metal ions with an 



orbital degeneracy are often relieved of this degeneracy 
by orbital ordering (due to Jahn- Teller distortions, for 
example). Orbital ordering can make neighboring ions 
inequivalent and often leads to the presence of both F 
and AF interactions in the nominally cubic perovskites. 
This is true for LaMnC>3, where a Jahn- Teller distortion 
elongates the oxygen octahedra and lifts the orbital de- 
generacy of the Mn 3+ t\ g e\ ion. Staggering of the elon- 
gated axis in the afe-plane minimizes strain and causes 
ordering of e g (3x 2 — r 2 /3y 2 — r 2 ) orbitals. In the ab- 
plane, e g orbitals on all neighboring ions have a half- 
filled/empty configuration, leading to F exchange, while 
overlaps of the ti g orbitals along the c-axis remain AF 
(half- filled/half- filled). The net result is the ^4-type mag- 
netic structure of ferromagnetic afe-planes coupled an- 
tiferromagnetically to neighboring planes along c. In 
LaVC>3, the degeneracy of the V 3+ t\ g ionic configura- 
tion is lifted by orbital ordering resulting in full occu- 
pancy of the xy-orbital and staggered occupancy of the 
xz/yz-orbitals in all three cubic directions. This orbital 
ordering leads to the C-type magnetic structure of AF 
planes coupled ferromagnetically along c. 

The energy of the exchange interactions in these and 
other perovskite magnets are most commonly determined 
by measurement of the spin wave dispersions using inelas- 
tic neutron scattering on single-crystal samples. Such 
measurements can be time consuming (taking a week, or 
more, of measurement time) and depend on the avail- 
ability of large single-crystals. In this manuscript, we 
show that the different exchange interactions in the ab- 
plane and along the c-axis (termed the magnetic ex- 
change anisotropy) in simple systems can be obtained 
rapidly by measurement of inelastic neutron scattering 
from powders. Inelastic neutron scattering from powder 
samples gives information on the spin wave spectrum re- 
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lated to the spin wave density-of-states (SWDOS). Such 
measurements can take as little as a few hours to com- 
plete and allow for efficient systematic studies of the de- 
pendence of exchange interactions on temperature and 
composition. 

We present powder inelastic neutron data for the afore- 
mentioned compounds; LaFeC>3, LaVC>3, and LaMnC>3, 
as representative of G-, C-, and A- type antiferromagnets, 
respectively. The data are compared to calculations of 
the spin wave spectra and their neutron scattering cross- 
sections using a Heisenberg model. The results show that 
the magnetic exchange anisotropy of F and AF inter- 
actions can be determined in a straightforward manner 
from powder data. 



II. SPIN WAVES IN CUBIC PEROVSKITES 

Spin waves in the cubic perovskite insulators can be 
described by the Heisenberg Hamiltonian. In the case 
where cubic symmetry is broken by charge or orbital or- 
derings, exchange interactions can become anisotropic. 
For the simplest kind of anisotropy, the exchange within 
the perovskite afo-plane (J a b) and that along the c-axis 
(J c ) have different values, and can even have different 
signs. The Heisenberg Hamiltonian becomes, 

H = —J a b J]] Si-Sj — J c Sj-S j—g^BH a s y^ajSi 

<i,j>\\a,b <i,j>\\c i 

(1) 

where S, is the spin vector on the ith site. The subscripts 
< i, j > || a, b (c) indicate that sums are restricted to near- 
est neighbor spins in the ab plane (along the c axis) . are 
restricted to nearest neighbor spins in the a6-plane and 
along the c-axis. Exchange energies are defined such that 
a positive J represents ferromagnetic exchange. Uniax- 
ial single-ion anisotropy is represented by an anisotropy 
field H a that acts on spin Si and points along the local 
spin direction (given by cti = ±1). 

In the following, we identify four different magnetic 
structures, one of which is ferromagnetic and the other 
three are antiferromagnetic varieties. The structures are 
differentiated by the signs of J a b and J c . 

F-type : J ab > 0, J c > 
G-type : J ab < 0, J c < 
C-type : J ab < 0, J c > 
A-type : J ab > 0, J c < (2) 

The spin wave dispersions for each type of magnetic 
ordering are obtained from a linear expansion (spin wave 
expansion) of the Heisenberg model4 When the single- 
ion anisotropy is zero, the dispersions are 



fiwip(q) = 25[2J o6 (l- 7+ (q))+J c (l-7,(q))] 
fiwc(q) - 2S{(2|J a6 | + |J c |) 2 

-[2|J ah | 7+ (q) + |J c | 7 ,(q)] 2 } 1 /2 

fiwc(q) = 2S{[2|J afc | + J c (1 - 7z (q))] 2 

-4J^(q)} 1/2 
hio A (q) = 2S{[2J ah (l- 7+ (q)) + |J c |] 2 

- Jhl (q)} 1/2 (3) 

where 7+ (q) = \ (cos q x a + cosq y a), 7 Z (q) = cos<j 2 a, 
q is the spin wave momentum, and a is the cubic per- 
ovskite lattice constant. The dispersions for each mag- 
netic structure are shown in fig. [1] for the case where 
\Jab\ — \Jc\- The notation for labeling the zone bound- 
ary reciprocal space positions are found in Kovalev.— The 
spin wave density-of-states (SWDOS) is the distribution 
of spin wave energies and is determined by the summa- 
tion over all wavevectors in the Brillouin zone (q), 

S H = 1^( W - W (q)) (4) 
q 

The densities-of-states are also shown for the four mag- 
netic structure types in fig. [TJ In addition, fig. [1] indi- 
cates the energies of the various extremal features in the 
SWDOS (van Hove singularities (vHs)) for any J ab and 

Jc- 
III. SPIN WAVE DENSITY-OF-STATES 

If \Jab\ = \Jc\i it is clear that the four types of ordering 
are easily discernible from the spin wave dispersions and 
densities-of-states. Figure [2] shows the SWDOS for each 
type of ordering in the case of | J ab \ — \Jc\- The max- 
imum spin wave energy increases as more ferromagnetic 
bonds are introduced. The F, A, C, and G-type struc- 
tures have the maximum spin wave energies of 12 J ab S, 
lOJabS, 8J a bS, and QJ ab S, respectively. 

In order to determine the exchange constants from 
powder samples, one must consider the degree of infor- 
mation available in the SWDOS. When \ J ab \ \ J C \, the 
positions of the vHs in the SWDOS allow identification 
of the exchange energies for each of the structure types: 

F-type: The SWDOS consists of five vHs as shown 
in fig. [Ha). The vHs corresponding to zone bound- 
ary spin waves along the c-axis or in the a6-plane lead 
directly to the corresponding exchange constants. Other 
zone boundary spin waves give sums of the exchange con- 
stants. FigureJ^a) shows how the spectral features of the 
F-type magnet develop as the ratio | J c /J a b\ changes. In 
the limit that J c = 0, the SWDOS of a 2-D ferromagnet 
consists of a sharp peak at 4| J a b\S and an upper cutoff 
at 8\J ab \S. 

A-type: The SWDOS consists of several vHs as indi- 
cated in fig. [Hb). The most useful for determining the 
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FIG. 2: (color online) Spin wave density-of-states for F-, A- 
C-, and G-type antiferromagnets in the case where | J a b \ = \Jc\ 
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FIG. 1: (color online) The spin wave dispersion along various 
symmetry directions (left panels) and the spin wave density 
of states (right panel) for (a) F-type, (b), A-type, (c), C-type, 
and (d) G-type perovskite antiferromagnets. Red arrows and 
labels indicate the energies of the extrema in the dispersion 
that give rise to van Hove singularities in the density-of-states. 



exchange constants (in units of \J a b \ S) are the maximum 
spinwave cutoff [2(4+ \ J C /J a b\)], the high energy edge of 
the central band [2(2 + \J C /J a b\)], and the low energy 
cusp (2 1 J c /Jab\)- Identification of any two of these vHs 
is sufficient to determine both J a b and J c . For example, 
J a b can be determined by the difference of the two high- 
est energy vHs identified here (equal to 2\J a b\S). Figure 
l^b) shows how the spectral features of the A-type anti- 
ferromagnet develop as the ratio \ J c j J a b\ changes. As J c 
becomes relatively small (in the approach to a 2D ferro- 



FIG. 3: (color online) Spin wave density-of-states for different 
ratios of the exchange \J a b/Jc\ for (a) F-type, (b), A-type, (c) 
C-type, and (d) G-type antiferromagnets. 

magnet), the central band collapses to a single energy at 
4| JaftlS 1 , the cusp just below the cutoff energy disappears, 
and the lowest energy cusp at 2 J C S disappears. This is of 
course identical to the same limit in the F-type magnet 
discussed above. 

C-type: The SWDOS consists of three vHs, as shown 
in fig. [IJc). The lowest energy peak is due to zone 
boundary spin waves in the basal plane and depends only 
on J a bS, while the splitting of the two main peaks in 
the SWDOS gives J C S directly. The development of 
the spectral features of the C-type antiferromagnet with 
the ratio |J c /J a t| are shown in figEJc). As J c becomes 
relatively small in the limit of the 2-D antiferromagnet, 



the two main peaks merge into a single peak at 4J ao S. 

G-type: The SWDOS consists of a single sharp peak 
at the cutoff energy with a very weak cusp-like vHs just 
below the cutoff energy. The energy of the peak in 
the DOS is determined by the average exchange (J) = 
(4 J Q 6 + 2J C ) /6. As shown in fig. [31(d) , varying the ra- 
tio J c / Jab shifts the entire spectrum. The weak second 
vHs is often masked by the finite energy resolution of the 
neutron spectrometer and in this case it is possible to 
determine only the average exchange < J >. 



IV. MEASUREMENTS 

In order to test the predictions made in the Heisen- 
berg model calculations above, inelastic neutron scatter- 
ing (INS) measurements were performed on the Pharos 
spectrometer at the Lujan Center of Los Alamos National 
Laboratory. Pharos is a direct geometry time-of-flight 
spectrometer and measures the scattered intensity over a 
wide range of energy transfers (hio) and angles between 
1° - 140° allowing determination of a large swath of the 
scattered intensity, S(Q,u>), as a function of momentum 
transfer (HQ) and fku. 

Powder samples of LaMn0 3 (LMO), LaV0 3 (LVO), 
and LaFeOs (LFO) were prepared by conventional solid- 
state reaction method and subsequently annealed to tune 
oxygen stoichiometry. Samples weighed approximately 
50 grams each and were characterized for phase purity 
by x-ray powder diffraction. Powders were packed in flat 
aluminum cans oriented at 45° or 135° to the incident 
neutron beam and INS spectra for LMO, LVO, and LFO, 
were measured with incident energies (Ei) of 75, 75, and 
160 meV, respectively. The time-of-flight data were re- 
duced into frw and scattering angle (29) histograms and 
corrections for detector efficiencies, empty can scattering, 
and instrumental background were performed. 

Unpolarized inelastic neutron scattering contains con- 
tributions from both magnetic and phonon scattering. 
In order to isolate the spin wave spectrum, the magnetic 
scattering must be separated from the phonon scatter- 
ing. This is accomplished by using the fact that the 
magnetic scattering falls off with Q (or 20) due to the 
magnetic form factor, while phonon scattering increases 
like Q 2 . Figure BJa) shows the full spectrum for LFO 
at T = 10 K as a function of 29 and Ku>, as reported 



TABLE I: Experimental conditions and data analysis pa- 
rameters for neutron scattering measurements on LaMn03, 
LaVC>3, and LaFe03 



Compound magnetic 


Ei (me 


V) low angle 


high angle 


ordering 




range 


range 


LaMnOa A-type 


75 


12-42° 


60-120° 


LaVC-3 C-type 


75 


7-32° 


60-110° 


LaFeOs G-type 


160 


1-31° 


55-95° 




-25 25 50 75 100 125 
Energy Transfer (meV) 



FIG. 4: (color online) (a) Inelastic neutron scattering inten- 
sity of LaFeOs (color scale) versus scattering angle and energy 
transfer at T = 10 K and Ei — 160 meV. Horizontal white 
lines delineate regions where phonon and magnetic scattering 
are isolated, (b) Neutron intensity summed over the angle 
range from 55 - 95° originating from phonons. (c) Neutron 
intensity summed over the low angle range from 1 - 30° (dots) 
and phonon background from scaled from high angle sum (ma- 
genta hatched region) (d) Isolated magnetic scattering from 
LFO (green) at T = 10 K. 



previously^ The band at 75 meV has intensity that falls 
off with 29, indicating that it is magnetic in origin. Data 
summed over the high angle range contains only phonon 
scattering (fig. 1Kb)), while the low angle range contains 
scattering from both phonons and spin waves (fig.[4jc)). 
The magnetic scattering component can be obtained by 
subtracting the high angle data from low angle data after 
scaling by a constant factor, as shown in fig. (He) for LFO. 
Figure BJd) shows that the resulting magnetic intensity 
for LFO indeed consists of a single peak at Hoj ~ 75 meV 
consistent with the G-type SWDOS shown in fig. [21(d). 

The strong peak at meV is elastic scattering and 
very weak peaks at ~20 and 30 meV arise from imper- 
fect phonon subtraction. The subtraction of the phonon 
intensity is subjected to error, primarily due to the fact 
that the scaling of phonon intensity from high angles to 
low angles is only expected to work for incoherent scatter- 
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FIG. 5: (color online) Extracted angle- averaged magnetic in- 
tensity (dots) versus the intensity calculated from a Heisen- 
berg model for the spin waves for (a) LaMn03, (b) LaV03, 
and (c) LaFeOs. 



V. CALCULATIONS OF THE SCATTERED 
INTENSITY 



When performing an INS experiment on a powder, the 
resulting INS intensities arise from the averaging of the 
inelastic scattering structure factor S(Q,u>) over all ori- 
entations of the crystallites. Despite the orientational 
averaging, the spectra can show evidence of the spin wave 
dispersions, especially at low angles (within the first Bril- 
louin zone) and in the vicinity of the first few magnetic 
Bragg peaks. Such dispersive features are clearly seen 
in the intensity plots of S(Q,ui) in figs [Ha), (c), and 
(e) for LMO, LVO, and LFO, respectively. Due to the 
weighting of the spin wave modes by coherent scatter- 
ing intensities, the Q-averaged intensity, S(to), as shown 
in fig. [5] does not necessarily give the SWDOS. This 
is only true in the incoherent scattering approximation, 
which does not apply to the case of scattering from a 
magnetically ordered system. Therefore, model calcu- 
lations of the powder averaged spin wave intensities are 
necessary for accurate comparison to the data. 

Numerical calculations of the spin waves in the linear 
approximation to the Heisenberg model give not only the 
dispersion relation cu n (q) for the n th (degenerate) branch 
(as shown in eqn. |4|), but also the spin wave eigenvec- 
tors, T„i(q), for the i th spin in the magnetic unit cell. 
The dispersion and associated eigenvectors can be used 
to calculate the spin wave structure factor for unpolar- 
ized neutron energy loss scattering from a single-crystal 
sample, S mas (Q,u/). 



WQ,c) = i( 7 r ) 2 (l + ^-^ 



E 



^^(Q)a 4 y^T m (q)e- 4Q - d ' 



(5) 



x (n(ui) + 1)6(lu - Lj n (q)) 



ing from a monatomic sample. For real multicomponent 
samples, the phonon intensity may not scale uniformly to 
low- Q due to coherent scattering effects (dependence of 
the phonon cross-section on Q) and also due to the differ- 
ent Debye- Waller factors for each component. It is diffi- 
cult to quantify this error without detailed phonon mod- 
els, however, based on the general agreement between 
the isolated magnetic scattering and the calculations dis- 
cussed below, the introduced errors are often small. 

The LMO and LVO data were treated in a similar- 
fashion to the LFO data. Figure [5] shows the isolated 
magnetic intensity for the three different antiferromag- 
nets. In each case, the magnetic spectra share similar 
features to the respective calculated SWDOS shown in 
figs. H and H 



where the i spin with magnitude Si pointed in direc- 
tion p, is located at position d^. crj = ±1 is the di- 
rection of the spin relative to the quantization axis /t 
for a collinear spin structure, q = Q — r is the spin 
wave wavevector in the first Brillouin zone. Finally, the 
function n (to) is the temperature dependent Bose factor 
and .Fj(Q) = ^gifi(Q)e~ Wi ^ is a product of the Lande 
g-factor, magnetic form factor, and Debye- Waller factor 
for the i th spin, respectively. The constant (^r ) = 
290.6 millibarns allows calculations of the cross-section to 
be reported in absolute units of (millibarns Steradian -1 
meV -1 (formula unit) For the simple perovskite mag- 
nets studied here, all ions in the magnetic cell are con- 
sidered to be equivalent. The structure factor can then 
be written, 
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TABLE II: Values of the ionic spin, exchange energies, and 
uniaxial anisotropy energies used in calculations of the neu- 
tron scattering intensity from spin waves. 



LaMnOo (A-type) 



Compound 


S 


Jab (meV) J c (meV) g^BH a (meV) 


LaMnOa 


2 


1.85 -1.1 0.6 


LaVOa 


1 


-7.8 2.9 0.6 


LaFeQ 3 


5/2 


-4.87 -4.87 



E 



^a l T m (q)e- l « d ' 



(6) 



x {n{uj) + \)8{uj - ct/„(q)) 



In the calculations below, we use the isotropic magnetic 
form factors found in the International Crystallography 
Tables^ and the Debye- Waller factor is set equal to one. 
The differential magnetic cross-section is proportional to 
the structure factor. 

To compare Heisenberg model spin wave results to the 
powder INS data, powder-averaging of S m ag{Q, is per- 
formed by Monte-Carlo integration over a large number 
of Q-vectors lying on a constant-Q sphere, giving the ori- 
cntationally averaged S mag (Q,uj) which depends only on 
the magnitude of Q. Figs [^b),(d), and (f) show calcu- 
lations of SmagiQ,^) (broadened by instrumental reso- 
lution) at 10 K for LMO, LVO, and LFO, respectively, 
and can be compared to the corresponding data in figs. 
Ilja), (c), and (e). In order to properly calculate the 
structure factor, we use the structural parameters for 
the three compounds (which are orthorhombically dis- 
torted perovskite structures). However, only the cubic 
exchange interactions J a b and J c are employed. The 
values for J a (,, J c , and H a used in the calculations for 
each compound are shown in Table [II] In the case of 
LMO, the anisotropy field was determined by Hirota et 
al. from single-crystal dispersion measurements and 
the exchange constants can be compared to the values 
obtained in that paper. For LVO, the anisotropy field 
was determined from cold neutron measurements of the 
anisotropy gap in powder samples^ The small anistropy 
energies reported here have very little effect on the energy 
of the zone boundary spin waves, which are determined 
primarily by the exchange. 

The calculations can be summed over scattering angles 
in order to compare the equivalent angle-summed data, 
as shown in fig [5] Overall, the agreement between the 
data and calculations is excellent. This is a testament 
to the effectiveness of the Heisenberg model for these 
compounds in predicting not only the spin wave ener- 
gies, but also the intensities. However, some differences 
observed in the comparison of data and calculation are 
worth noting. At low energies near to the elastic line, 
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Momentum Transfer (A" 1 ) 

FIG. 6: Panels show the measured (left) and calculated (right) 
neutron intensities for (a)-(b) A-type LaMn03, (c)-(d) C-type 
LaVOs, and (e)-(f) G-type LaFe03. Experimental conditions 
and calculation parameters are given in the text. For the 
measured data in panels (a), (c), and (e), phonon intensities 
have not been subtracted and lead to a more complicated 
spectral image as discussed in the text. For each panel, the 
curved white lines indicate the low-angle summation regions 
leading to the corrected magnetic spectra in fig. 5. 



additional intensity is observed, most notably in the 5 
- 10 meV range in Ei =75 meV data. The origin of 
this intensity is unclear, but it is possible that it origi- 
nates from multiple elastic scattering. For LaMnOa, this 
additional scattered intensity, combined with insufficient 
elastic energy resolution, does not allow the observation 
of the low energy vHs expected at ~6 meV that can be 
used to determine J c . Higher resolution measurements 
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are required to obtain J c exclusively. 

The poorest agreement between data and calculation 
occurs for LaVC>3 in fig^b). While the data shows clear 
vHs at ^32 and ~44 meV, the calculated intensity shows 
a shoulder, rather than a clear peak in the upper vHs, 
indicating that the Heisenberg model does not reproduce 
the LVO spin wave intensities with the same accuracy as 
for LMO and LFO. Much of this discrepancy may be due 
to the physics of LVO, where C-type magnetic strucure 
can arise from either weak Jahn- Teller driven orbital or- 
dering or orbital singlet formation. Competition between 
these two scenarios require spin-orbital coupling terms 
that go beyond the Heisenberg Hamiltoniani 10 ' 11 Also, 
additional intensity at 20 meV and 50 meV in LVO is 
likely due to improper phonon subtraction. This sub- 
traction is made more difficult due to the small spin 
(S — 1) of the V ion, which leads to weaker spin wave 
scattering (see eqn. [6]). 

Figure [5] shows a comparison of the measured and cal- 
culated scattered intensities as a function of Q and fku 
for the three compounds. The calculation results in figs. 
Eth), (d), and (f) show clearly the coherent scattering 
of the powder averaged spin waves. The most obvious 
coherent scattering feature is the necking down of acous- 
tic spin waves in the vicinity of allowed magnetic Bragg 
reflections. The characteristic ordering wavevectors for 
the different antiferromagnets are; (0,0,1/2) for A-type, 
(1/2,1/2,0) for C-type, and (1/2,1/2,1/2) for G-type (us- 
ing the cubic indexing) . The first two observed magnetic 
Bragg peaks in each case are indicated in fig. [5J Ad- 
ditional coherent scattering features can also be seen for 
zone boundary spin waves, where intensities tend to peak 
in between the allowed magnetic Bragg peaks. Fig. [6] 
enforces the general agreement of the Heisenberg model 
calculations of the spin wave intensity with neutron scat- 
tering measurements. Unfortunately, the comparison of 
the Heisenberg model spin wave intensities to the data 
is complicated because measurements also contain coher- 
ent phonon scattering intensity. The phonon intensity 
bands present themselves mainly as horizontal (constant 
energy) streaks. A prominent phonon band can be seen, 
for example, at 25 meV in LVO (figH^c)), and at 25, 40, 
and 60 meV in LFO ((fig He)). 

The success of the Heisenberg model in estimating 
the measured spin wave intensities is better observed 
by plotting constant energy Q-cuts, as shown in fig. 
for LMO. The plots show Q-oscillations of the exper- 
imental magnetic spin wave scattering above a back- 
ground due mainly to phonon scattering and back- 
ground/multiple scattering. A constant background and 
incoherent phonon scattering intensity (proportional to 
Q 2 ) are added to the calculated spin wave scattering in 
order to compare to the measured data. The agree- 
ment is excellent. The overall consensus is that the 
spin wave intensities arc well represented by the Heisen- 
berg model and it is promising that one can obtain more 
from powder data than just an estimate of the spin wave 
DOS. Analysis of the full structure factor S mag (Q,u>) 
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FIG. 7: (color online) The Q-dependence of the neutron scat- 
tering data for different energy transfer ranges in LaMnOg: 
(a) 8-12 meV, (b) 12-16 meV, (c) 18-22 meV, and (d) 26-30 
meV. The red dots are the experimental data. The dashed 
lines are estimates of the incoherent phonon background plus 
multiple scattering. Solid lines are calculations of the poly- 
crystalline averaged spin wave scattering using the parameters 
in the text plus background. 



may allow exchange interactions to be determined in 
more complicated magnetic structures, or with interac- 
tions beyond nearest-neighbor. However, a full analy- 
sis of powder averaged spin waves requires better un- 
derstanding of the phonon spectra and multiple scatter- 
ing. In the future, we plan on combining fully coher- 
ent calculations of both phonons and spin waves to at- 
tempt a more ambitious analysis of the full S{Q,uj) = 



Smag(Q, + S p h 



12 



VI. SUMMARY 

We have demonstrated that inelastic neutron scatter- 
ing experiments on powders, in combination with calcu- 
lations of the spin wave scattering in a Heisenberg model 
can give detailed information about the exchange inter- 
actions in simple magnets. The agreement of not only 
the spin wave DOS, but also the Q-dependence of coher- 
ent features in the spin wave scattering gives hope that 
even more complicated magnetic systems can be analyzed 
using the full S(Q,lu). The advent of new spallation 
neutron sources, such as the Spallation Neutron Source, 
will allow the rapid measurements of samples with good 
statistics and make detailed systematic studies of mag- 
netism possible. 
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